CERN-TH/2002-195 
UTHEP-02-0802 
MPI-PhT-2002-036 



Electric Charge Screening Effect in 
Single-W Production with the KoralW Monte Carlo* 

S. Jadach a , W. Placzek fec , M. Skrzypek ac , B.F.L. Ward dec and Z. Was a c 

"■Institute of Nuclear Physics, ul. Radzikowskiego 152, 31-342 Cracow, Poland 
b Institute of Computer Science, Jagellonian University, 
ul. Nawojki 11, Cracow, Poland 
C CERN, Theory Division, CH-1211 Geneva 23, Switzerland 
d Department of Physics and Astronomy, 
The University of Tennessee, Knoxville, TN 37996-1200, USA 
e Werner-Heisenberg-Institut, Max- Planck- Institut fur Physik, Fdhringer Ring 6, 

80805 Munich, Germany 

Abstract 

Any Monte Carlo event generator in which only initial state radiation (ISR) is 
implemented, or ISR is simulated independently of the final state radiation (FSR), 
may feature too many photons with large transverse momenta, which deform the 
topology of events and result in too strong an overall energy loss due to ISR. This 
overproduction of ISR photons happens in the presence of the final state particle 
close to the beam particle of the same electric charge. It is often said that the lack 
of the electric charge screening effect between ISR and FSR is responsible for the 
above pathology in ISR. We present an elegant approximate method of curing the 
above problem, without actually reinstalling FSR. The method provides theoretical 
predictions of modest precision: < 2%. It is, however, sufficient for the current 1W 
data analysis at the LEP2 collider. Contrary to alternative methods implemented 
in other MC programs, our method works for the ISR multiphotons with finite pt- 
Although this method is not an exact implementation of the complete/exact ISR, 
FSR and their interference, it is very closely modelled on it. We present a variety 
of numerical results obtained with the newest version of the KoralW Monte Carlo, 
in which this method is already implemented. 
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1 Introduction 



In the Monte Carlo programs for e~e + colliders, one has to model initial state radia- 
tion (ISR) as precisely as possible. For some processes where the precision requirements 
are moderate, such as the single-W^ production process at LEP, where the combined 
LEP2 precision is estimated at ~ 7%, see Ref. |l|], it is sufficient to implement ISR in 
the leading-logarithmic (LL) approximation. This can be done by modelling ISR in the 
collinear, p T = 0, approximation (see for instance refs. @,[§), by "unfolding" the strictly 
collinear structure functions 0] or using the so-called "parton shower" technique as in 
Ref. || (recently employed also in 0), or the Yennie-Frautschi-Suura (YFS) exclusive 
exponentiation employed in Ref. f?j]; see Ref. [f| for a more complete list of references. In 
all of the above MC programs the final state radiation (FSR) from charged final particles 
is either not modelled at all or simulated completely independently of the ISR, which is 
the general spirit of the LL approximation. This may be problematic if one of the final 
state charged particles gets close to the beam particle with the same electric charge. In 
such a case ISR photon emission is damped, because of the well-known electric charge 
screening effect (ECS), and the total energy distribution of the ISR photons gets softer. In 
the terminology of the LL it is described as the "decrease of the LL scale" from s, which is 
the square of the centre-of-mass-system (CMS) energy, to t, which is the four-momentum 
transfer (squared) from the beam to the nearby particle of the same charge. We shall 
refer to the above phenomenon as the leading-logarithmic scale transmutation (LLST). 

Let us stress that in the calculations based from the beginning on the YFS exclu- 
sive exponentiation, as in the BHWIDE or KKMC g§|TI| MC programs, the ECS 
and LLST are automatically built in to an infinite order, hence there is no need to re- 
install them. The ECS and LLST in exclusive exponentiation are direct manifestations 
of the interference effects between ISR and FSR. Hence, they are also often described as 
"coherence effects". 

As already said, owing to the lack of ECS in the Monte Carlo calculation in which 
ISR is modelled using the s-scale, one may encounter two pathologies: (a) overproduction 
of photons with large transverse momenta, which results in the noticeable deformation of 
the topology of events and (b) too strong an overall energy loss due to ISR. In particular 
this happens in the single-W 7 production process at LEP2, when modelled by the 



MC program KoralW p2|, |13|, [L4fl . The aim of this work is to invent a simple method 
of correcting for ECS, and curing the above deficiencies, without the need of generating 
FSR in the MCQ. The method is quite general, although its most immediate application is 
improving the theoretical prediction for the 1W production process at LEP2 to the < 2% 
precision level. Contrary to alternative curing methods implemented in the other MC 
programs P,f|, our method works for finite-pr ISR multiphoton distributions. It is not 
an "ad hoc" method because its development is based on the consideration of the exact 



1 An obviously better alternative is to implement YFS exclusive exponentiation for the entire e~e + — > 
4/ process. This is however too much work to be completed before the end of the LEP2 data analysis. 
For the e~e+ — > WW —> 4/ subprocess, where precision requirements are ~ 0.5%, the ECS/coherence 
effects between ISR and FSR (W decays) are strongly suppressed. 
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YFS exclusive exponentiation for ISR, FSR and their interference |TTJ . 

After defining our method, we shall also present a variety of numerical results exploit- 
ing the newest version of the KoralW Monte Carlo, in which this method of introducing 
ECS is already implemented. 

The electric charge screening effect together with leading-logarithmic scale transmu- 
tation for the multiphoton emission are described in the following section. The proposed 
method is practical because: 

• the effects are introduced by means of a simple well-behaving Monte Carlo weight, 

• it does not require generation of the additional photons in the final state (FSR), 

• it provides an ISR precision good enough for LEP2 data analysis, that is < 2%. 



The construction of the ECS weight is first described in Section |1.2| , using a simple 
example of the process — > for which the ECS effect and the method of 

its introduction in the ISR is illustrated with the help of many examples of the photon 
distributions. This method is later generalized to the e~e + — > 4/ process in Section [| 
where also the relevant correction to the virtual form factor (overall normalization) is 
defined. Numerical cross-checks and final discussion are provided in Section 01. 



1.1 Electric charge screening 

The electric charge screening effect is a well-known phenomenon of suppression of the 
photon emission from any subset of compensating charges (charges adding to zero or to 
small total charge), which are close in the momentum phase space, i.e. which have small 
effective mass, small angular distance, etc. One should exclude from the consideration 
charges/particles that are well time- separated, i.e. from the decay and production parts 
of the narrow resonance production/decay processes. A good example of a subset of 
compensating charges is for instance a fast-moving light fi + fi~ pair - a "dipole" of com- 
pensating charges. For the compensation considerations, electric charges of every initial 
state particle should be assigned an additional minus sign with respect to particle charges 
in the final state (much as in the definition of the four-momentum transfer). For instance 
in the low-angle Bhabha (LABH) process, the initial state e~ and the final state ej form 
together a "dipole" of compensating charges and the ECS/LLST takes place, which means 
that the emission of photons with an angle greater than the angular size of the (e~,ej) 
compensating system gets strongly suppressed. 

The leading-logarithmic scale transmutation is just another face of the same ECS. For 
ISR taken alone, the LL parameter L e (s) = 2^1n^|, which controls the strength of the 

loss of the total energy due to ISR, is defined for the so-called LL energy scale y/s, being 
just the total centre-of-mass system (CMS) energy. Actually this so-called big logarithm 
L e (s) comes from the integral over the photon angle L e (s) ~ J^ 2 , Taking into account 
the interference between ISR and FSR introduces ECS. This becomes important when, 
for instance, one of the final state e~'s approaches the initial beam e~ by an angle 9. 
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Then, the photon angular distribution gets suppressed (i.e. loses the l/i? dependence) for 
d > 6. Hence L e (s) gets replaced in the so-called LL ISR structure functions by a smaller 
L e (t) ~ fjrf/s IF- Obviously, LLST is a numerically important effect when \t\ <C s. Such 
is the case of the 1W production, where at least one of the final state electrons gets inside 
the beam pipe due to the strong peak in its scattering angle (because of t-channel photon 
exchange in the Born distribution). 

All the above discussion was for real photons, and in fact the whole ECS can be 
viewed as an almost completely classical phenomenon. The real photon angular distri- 
bution, in the soft-limit approximation, is just the square of the classical electric current 
vector. Quantum-mechanical virtual corrections merely correct the overall normalization, 
cancelling properly infrared (IR) singularities in the photon energy distribution. (The 
IR-divergent part of the exponentiated virtual correction has necessarily L e (t) as a coef- 
ficient.) 



1.2 "Mhamha" toy example 

As already indicated, in our method of implementing the ECS effect in ISR we do not 
include explicit FSR; however, we take as a guide the YFS model in which ISR, FSR and 
their interference are included. In this section we examine carefully the ECS effect for a 
single real photon in such a complete model, proposing at the end of the exercise a simple 
correcting weight introducing the ECS effect in the pure ISR case. 

The above will be done by carefully scrutinizing several variants of the angular photon 
distributions of a "toy process" : 

H~(Pa) + V + (Pb) -> ^~(Pc) + V + (Pd) + j(k), (1) 

a muonic analogue of the Bhabha process, which we shall call a "Mhamha scattering 
process"^. In the soft photon limit, the final state distribution is equal to the Born 
distribution times the complete soft photon factor 

§ bd = a f Pa Pb Pc | Pd V ^ 

47T 2 \kp a kp b kp c kp d J 



the same factor which is a basic element in the YFS exponentiation [15]. In order to see 
the ECS effect we compare the above complete S a bcd with the incoherent sum S a b -\~ Sd 
of the ISR and FSR contributions, neglecting the interference IFI=ISR®FSR, where 

g = a_ (Pa_ _ Pb_\ 2 g = a_ (Pc_ _ Pd_\ 2 

47r 2 \kp a kp b J ' cd Att 2 \kp c kp d J 

We choose a relatively small scattering angle 9 a b- In Fig. [l]we plot S a b + S c d and complete 
Sabcd side by side, as a function of the photon polar variables (cosi?, <p). In order to see very 
clearly the structure of the distribution over the entire unit sphere, we choose (cos $,</?) 

2 The name "Mhamha" was suggested by A. Blondel in the context of the muon collider studies. 
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with respect to the z-axis pointing perpendicularly to the beam axis (beams are along 
the x-axis). We choose \/s = 5 GeV at which muons are already very relativistic and the 
muon scattering angle 20° - small enough to see ECS effect. The four peaks of the photon 
emission intensity in Fig. [1] are centred around directions of the four charged muons. The 
ECS effect is seen as a clear suppression of the photon emission intensity beyond the two 
dipole peaks (a, c) and (b,d). The first dipole is for the incoming and outgoing \T and 
the second one is for the incoming and outgoing fi + . In Fig. [I] the two t-channel dipole 
peaks are much sharper for the complete S a b c d than for the incoherent superposition of 
ISR and FSR, S a b + S c d- The top part of the dipole peaks is worth inspection, see Fig. 0, 
as it features the well-known "helicity zeros" as deep "craters", exactly in the direction 
of the emitter charge. The craters are of an angular size m^/Ebeam and will shrink to 
a negligible size at very high energies. The ECS effect is also very clearly seen again in 
Fig. |2], as a strong sharpening of the photon radiation intensity beyond the dipole double 
peak. 

Is the radiation pattern of S a bcd, which we see in Fig. specific to a t-channel character 
of the process? Not really. The radiation pattern of the real /J,~fi + pair with effective mass 
equal to y/\t\ and boosted to the total energy of 5 GeV is completely undistinguishable 
from what we see in Fig. || The other way of describing it is that S a bcd is extremely well 
approximated by the incoherent sum of S ac + Sbd- Let us stress that S ac and Sbd look 
very simple and natural in the two Breit frames (rest frames of (a, c) and (b, d)) and the 
non-trivial shape of the photon intensity in Fig. |^ is the result of the trivial Lorentz boost 
from the Breit to the CMS frame. (In a sense, the non-trivial shape plotted in Fig. |2| is 
the "fault" of the observer himself, who has chosen to examine it from an "unnatural" 
reference frame.) The other interesting observation is that if we have plotted S ac + Sbd 
in Figs. [I] and [| then we will have found that it is completely undistinguishable from the 
Sa&cdQ- In order to see any noticeable difference between these two distributions, we have 
to go to the backward scattering angle. We do this in Fig. ^|, where we compare these 
two emission distributions for the scattering angle 9 = 135°. The peaking structure is the 
same; however, for S a b c d we see a rich interference pattern, especially for photons far away 
from the four charges. At this large scattering angle we have checked that the S ac + Sbd 
and S a b + Scd distributions look almost identical. 



1.3 ECS weight — preliminaries 

Having seen in Figs. |I] and ^| what the coherence ECS effect is, let us now ask the important 
practical question: Provided that, as in KoralW, we have only ISR, could we modify the 
ISR distribution in such a way that we get the same sharpening of the photon angular 
distribution, as in the real world of the coherent S a baP- Our proposal for the correcting 
weight is the following: 

*W*) = , (4) 

Sab(k) + S cd {k) 



3 This is why the BHWIDE MC program, which models S ac + Sbd at the low MC level, is so efficient 
in modelling the multiphoton radiation in the complete wide-angle Bhabha process. 
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and we are going to show that it does what we want. Let us examine the corrected ISR 
distribution 

S ab (k)W EC s(k). (5) 

First, we notice that W^csik) — 1 (apart from the unimportant helicity zero) close to the 
angular position of the outgoing particles, so it does not try to re-install the FSR, and the 
corresponding MC weight should not have an inconvenient long tail. Then, what is most 
important, Wecs really cuts off very strongly photon emission at angles greater than the 
scattering angle 8 ac or 9m- This is seen very clearly in Fig. |], where the "angular size" of 
the corrected ISR emission pattern shrinks by a factor of ~ 3 with respect to the original 
one. This is clearly visible even for our moderately small scattering angle 9 ac = 20°. In 
the real KoralW application we shall, of course, apply the product of the Hecs(&) over 
all ISR photons. The overall normalization will also be corrected by means of correcting 
the exponential virtual+real form factor, which also compensates for the IR divergence 
of the average ECS weight dlj ^ / ecs(^))- 

How does one justify the fact that we apply a correcting weight, which includes 
ISR®FSR interference, without generating FSR in the MC? A quite general answer is 
that the above can be understood as a procedure of integrating analytically over the real 
FSR photons and including the result of the integration in the normalization^. Such a 
procedure can be very close to experimental reality, provided the FSR photons are com- 
bined with the outgoing charged particles. From the theoretical point of view the above 
procedure is quite common in the world of the LL approximation. For example, a similar 
solution is used to install interference between different branches of FSR in the PHOTOS 
Monte Carlo As we shall see later on, our procedure not only makes sense in the 



LL approximation, but also features the precise IR cancellations and has the expected LL 
behaviour. 



2 ECS Correction weight for the 4/ process 

We have already explained what the ECS effect is, and sketched how to introduce it in an 
approximate way in the ISR Monte Carlo, with the help of the correcting weight Wecs- 
All that was done for the 2f final state in a rather qualitative way. In the following we 
define the analogous correcting procedure for the four-fermion (4/) final state; we shall 
also present several quantitative tests, see Section [|. 

2.1 Real emission part of Wecs 

Let us start by defining the real emission part of the W^ecs f° r the 4/ process e~e + — > 
fc{p c ) + IdiPd) + feiPe) +ff(Pf)- F° r this process the complete soft photon emission factor 

4 The FSR photons can be added later on; however, care must be taken to get the ECS for the FSR 
in a similar way as for ISR. 
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reads 

^ = lE J^, = -j^)', (6) 

where Z^ = ZiZjOiOj, Zi is the sign of the i-th charge, and 9i = + (— ) if particle i is 
outgoing (incoming). The natural extension of the weight of Eq. (|j) would be 

Sab(k) + Scdefik) 1 

where S c d e f is defined as in Eq. (H), restricting summation to the FSR part, i,j = c, cf, e, /. 
The above weight could be implemented in the MC without much problem; however, 
it would complicate the construction of the accompanying normalization weight WQq™ 
described in the next section. 

We have noticed, however, that since we are working in the LL framework, we can 
simplify the weight of Eq. (|7|), and effectively replace it with the variant of the 2f weight 
of Eq. (f|). How is this possible? One has to keep in mind that our aim is to deal with the 
situation in which one of the final state particles (electrons) is close to the beam or two 
particles (electrons/positrons) are close to the beams. Let us call the above two situations 
(following a long established terminology) singly- and doubly-peripheral scattering. 

In the singly-peripheral (SP) situation we expect ECS in only one hemisphere. Photons 
will be emitted inside a narrowly collimated t-channel dipole close to the beam; there will 
be a gap in the photon angular distribution extending from the dipole down to the nearest 
particle in the "central region"^ out of the remaining three final particles. In the doubly- 
peripheral (DP) configuration, we shall have two narrowly collimated t-channel dipoles 
close to the beams and two gaps down to the central region. In fact the two gaps will join 
together in the central region. This gap structure will be clearly seen in the numerical 
results in Section |], where we shall plot the photon distribution in the rapidity variable 
y = -lntan(tf 7 /2). 

Now, if our main aim is to reproduce one or two of the ECS gaps in ISR radiation, 
and if we restrict ourselves to the LL approximation, then it does not matter how many 
particles there are in the central region. What really matters is how many (one for SP and 
two for DP case) particles close to the beams there are, and at which angles those closest 
to the beams (members of t-channel dipoles) are. Consequently, in the LL approximation, 
for the purpose of reproducing ECS, it is perfectly safe to replace the four final particles 
by just two (especially that there are only two beams). In the DP case, the choice is 
clear: we take those two particles closest to the beams and properly match the beam 
charges (form together with the beams t-channel dipoles), ignoring the other two (charge 
conservation is of course respected). In the SP case one charge is again the member of 
the t-channel dipole and another charge we place anywhere in the central regions; it may 
be tied up with one of the three remaining particles or not - it does not matter. 

5 By central region we understand angles greater than, say, 10° from both beams. 
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(7) 



The two final "effective final state particles" entering the ECS weight we shall denote 
C and D, and the real-emission part of our ECS weight we define as follows 



S ab {k) + S CD {k) + S aC (k) + S bD (k) + S aD (k) + S bC (k) 
ECS \ l 1 lh {) S ab (k) + S CD (k) 

(8) 

The above formula, apart from notation, is identical to that of Eq. (HI). Also, in the 



numerator we have used the algebraic identity of Ref. [15] in order to decompose S ab CD 
into six dipole-type contributions. 

In practice, the final state fermions C and D in Eq. ([8]) are determined by the kine- 
matics of the final state (FS). In the case of the so-called "single- W" (1W) final state, in 
which we are mainly interested here, we check whether one final state electron (positron) 
is lost in the beam pipe and the decay products of the 1W are visible outside the beam 
pipe. In such a case we identify the in-beam-pipe particle as one of the effective particles 
C or D (according to its charge). As for the other one, there is freedom of its choice, as 
explained above. Since all other FS fermions are at large angles (central region) it should 
not matter which one we pick. In fact we do something even more primitive: we assign 
the second remaining C or D as directed perpendicularly to the beam pipe. In the second, 
DP case of the e~e + ff final state with both e~e + lost in the beam pipe, we assign C and 
D to e~ and e + , of course. 

Let us summarize the main points on the weight of Eq. (H), as implemented in KoralW: 

• The only purpose of the weight W^cs * s ^° res tore the ECS effect due to ISR®FSR 
interference. 

• We do not aim at re-creating the FSR. This would be formally possible with a similar 
weight; however, it would lead to an awful weight distribution and a non- convergent 
MC calculation. 

• We get Wecs — ► 1 for photons collinear with the FS effective fermions C and D. 
This ensures a very good weight distribution. 

• The FSR can be treated separately, either inclusively (calorimetric acceptance) or 
exclusively, generated with the help of PHOTOSQ. 

Finally, we notice that in the ECS weight we may insert massless four-momenta with- 
out any problem - in practice we shall use a massless limit of Eq. (^): 



3 Care has to be taken to implement ECS for FSR, if necessary. 
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2.2 Virtual+soft correction to normalization 

The average of the real emission weight (W^cs) taken alone is infrared (IR)-divergent as 
~ In e, where e defines the infrared cut E 1 > eE^ am on the photon momentum in the CMS. 
In the YFS exponentiation of ISR in KoralW, based on Ref. |17| , the IR cancellations 
occur numerically between the real soft-photon factors Yl i S(ki) and the YFS form factor 
Fyfs = exp(2a£> + B(e)) (contrary to typical parton-shower MCs, where they are built- 
in features of the MC algorithm). Since it is well known how these IR cancellations do 



occur, see refs. [O, |TT|], it is therefore also possible to calculate the IR violation of the 
(^ecs)' an d to correct for the above divergence precisely. At the same time, we shall also 
correct for the LLST in the virtual part of the YFS form factor and the non-infrared LL 
corrections to the 0(a). As a result, the corrected total cross section will become again 
independent of the IR parameter e, as the original one. This we shall check numerically. 
Again, all of this procedure should be regarded as a "shortcut" with respect to the full 
YFS exponentiation with ISR+FSR including ISR®FSR interference, keeping in mind 
that our precision tag is limited to 2%, and that simplicity of the method is also our high 
priority. 

The correction to the overall normalization, which corrects for the IR divergence ex- 
actly, and introduces LLST in the virtual LL correction, reads as follows: 



W n E Zs = exp (j(7* - 7*)) exp (AU(e)) 



f d*k f d*k (10) 

AU(e)=U(e)-U R (e), U(e) = J j^S ab (k), U R (e) = J —S ab (k)w R (k). 

The factor exp(AtZ) cancels exactly the e-dependence and compensates approximately for 
the normalization change due to the (W^cs) weight. The proof and details of construction 
of exp (At/) are given in Appendix [A|. 

The factor exp((3/4)(7t — 7 S )) provides the LLST in the virtual part of the form 
factor. The 7 S is defined as 7^ = 2(a/7r)(log(s/mg) — 1). The scale of the 7* depends on 
the number of final state particles that effectively contribute to the interference weight 
(H). When only one line is modified, say e~, we have 7* = (1/2) (7s + 7t-), where 74- = 
2(a/7t)(log(\t~\/ml) — 1) and t~ denotes the square of the four-momentum transfer from 
the e~ line. When both e~ and e + lines are modified, we get j t = (1/2) (j t - + 7 t +). 

The triple integral of Eq. (TO) has to be computed for every generated event. It 
is therefore of high importance to be able to do it very fast, so that the Monte Carlo 
generation is not slowed down too much. We have not attempted to compute it completely 
in an analytic way. Instead, we perform two out of three integrations analytically and 
the remaining one numerically. Such a procedure proved to be fast enough. We show in 
Appendix [B| how to reduce the integral in Eq. (|T0| ) to a simple one- dimensional integral. 
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3 Running QED Coupling Constant 



In the practical applications, a more precise prediction for the lW^-type processes requires 
also the inclusion of the effect of the running QED coupling constant in the e~ (e + ) vertex 
from the value at scale to the actual small transfer value t~ (t + ). We include this 
effect in the form of an overall factor that multiplies the whole matrix element squared. It 
is activated only for the electron or positron line (in the Feynman diagram) for which the 
ECS corrective weight is actually applied. In the SP case of one e ± line, the correcting 
weight reads: 

'Run - W Run 



w Run = w± = I ] , fir 



where ac M is the value of coupling constant in the G M scheme^]. For the DP case (two 
lines) we put 

WRun = W~ Run W+ un . (12) 

The precision of such a naive solution is, in most cases (including leptonic final states), 
better than 2%, and for semileptonic final states even better than 1%, as discussed at 
length in Ref. W~ 



4 Numerical results 

In this section we check very carefully that our ECS weight 

W ECS = W n E Zs W r E e Ss W Run (13) 

introduces the ECS effects in the photon angular distribution, the LLST in the total 
energy loss due to ISR photons, and corrects the QED coupling constant for the hard 
process. 

4.1 General consistency tests 

In this subsection we present some numerical results that demonstrate the action of the 
weight emulating the ECS effect. 

The first test, necessary for consistency, shows that the cross-section is independent 
of the dummy IR cut-off e of Eq. QTO) . In Table 1 we show the values of the total cross 
sections for the ev e ud channel for a few values of the dummy IR cut-off e. The energy 
is set to 190 GeV, and the cut-off on the maximal angle of the scattered electron is 2.5° 
w.r.t. the electron beam in the effective CMS frame of outgoing particles. One can see 
that, within the statistical errors, there is indeed no dependence on the e. 

Having checked the self-consistency of the emulation, we can now proceed to demon- 
strate the main result of this work - the ECS suppression of the transverse radiation. In 

7 If any scheme other than G M is used in KoralW, the program will report a conflict and stop at this 
point. 
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Table 1: Value of the total cross section for the ev e ud final state, for various values of the dummy 
IR cut-off ejR. 



Fig. |5] we show (in the doubly-logarithmic scale) the differential distributions da/dy with 
(red dots) and without (blue open squares) the correction weight for the ECS effect for 
the e + e~ss channel. The angles of scattered electron and positron are set to be between 
2 and 0.2° w.r.t. the beam line (the double "ice-cream-cone" configuration around the 
beam line) in the laboratory frame. The distribution is obtained by summing over all 
photons of a given event and is normalized to arbitrary units. The y = — log 10 tan(# 7 /2) 
is proportional to the (pseudo)rapidity of the photon. 

The scattering cone of electron (positron) corresponds, in the y scale, to the values 
from 1.76 to 2.76 (—1.76 to —2.76) and are marked in the plot with blue shadow bands. In 
the Figure one can clearly see the ECS suppression between these bands. The spectrum 
for photon emission outside the area (i.e. inside the e~ and e + cones) remains flat and 
unchanged. Between the bands (i.e. outside the e~ and e + cones) the suppression is clearly 
visible. We have fitted the suppressed spectrum with straight lines. The values of the 
fitted slopes are shown in the figure to be in good agreement with the values of ±2, which 
means that we see the 1/0^ suppression of the photons beyond the dipole angular size, 
naively expected from the "multipole expansion" . 

In the next step we look into the IW-type final state of ev e ud. This time, there is 
only one electron in the final state, so we expect the suppression area to be asymmetric 
- only in the forward direction. The result of the simulation is shown in Fig. ^. In blue 
(open squares) we show the same curve as in Fig. |5], i.e. the e~e + ss case (appropriately 
renormalized) , whereas in red (dots) we present the new eu e ud result. Shadowed bands 
again visualize the location of the electron (positron) "ice-cream-cones". For the ev e ud 
case the area corresponds to the "ice-cream-cone" between 0.2 and 0.5°. The other end of 
the ev e ud suppression area is located at zero (90°), where the fictitious particle is located 
and leaves the radiation in the backward hemisphere unchanged. Note that the exact 
location of this fictitious particle (in the area of reasonably large angles) is irrelevant 
anyway, as the effect depends on it logarithmically. The fitted lines can be used as before 
to confirm the correctness of the suppression factor (slope). 

Finally, in Fig. [F] we show the ev e ud process for two different values of the half-opening 
angle of the electron "ice-cream-cones": 0.2 to 0.5° and 0.02 to 0.05°. One can see that the 
left (backward) edge of the suppression area is, as expected, identical in the two cases, 
whereas the right (forward) edge follows the electron opening angle (the range of this 
angle is marked with shadowed bands). 

Another characteristic distribution is the da/d\og 10 v distribution, with the definition 
v = 1 — s'/s. As this distribution is sensitive only to the value of the transfer in the 
hard scattering process, we have defined the acceptance cuts for this exercise to be —9 > 
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log(t/s) > —11. In the limit of small v (soft limit) we know from YFS exponentiation 
thatQ da/dv = exp(3/47)F(7)7^- 1 a((l - v)s) where F( 7 ) = exp(-7CW er )/r(l + 7) = 
1 + C(7 2 ). The 7 for the case of pure s-channel ISR is governed by the scale s, i.e. 7 S = 
2a/n(\og(s/ml) — 1). Inclusion of the ECS suppression leads to the replacement of the 
scale s — > \t\, i.e. 74 = 2a/ii(\og(\t\/ml) — 1). In the case of the lW-type process (ei> e ud) 
this suppression happens only in the electron line, whereas for the positron line we still 
retain the scale s, so the effective j t becomes % = 1/2(7* + It)- For the e + e~ case, the 
suppression happens for both lines and we have j t = It- 

In Fig. we show a ratio of da/dlog l0 v distributions with and without the ECS 
correction for the processes e + e~ss (blue open squares) and for ev e ud (red dots). For the 
e + e~ ss process we applied an additional cut-off of 3.5 mrad for the minimal angle of s and 
s quarks with respect to the beam line and we used the so-called "extrapolation procedure" 
that preserves the value of smallest transfer (see KoralW manual J?|] for details). As follows 
from the above discussion, we expect to see (in doubly-logarithmic scale) straight lines of 
the form a + h log 10 (w) where the slope of the line should hence be given by bth = Jt — Is = 
a/irlog(t/s) and the free coefficient by a th = 3/4(7 t - 7 S ) log 10 (e) + log 10 (7*/ 7 s)- m our 
exercise the average value of the \og(t/s) as calculated in the simulation is \og(to/s) = 
—9.86. This leads, at t = t , to the expected theoretical values of slopes of b t h = —0.023 
and a t h = —0.104 for tv e ud and b t h = —0.046 and a t h = —0.241 for e + e~ss. The 
corresponding results of the fits of the actual Monte Carlo simulations to straight lines, 
shown in the insets in Fig. ||, are b t h = —0.021 ±0.002 and a th = —0.093 ±0.005 for tv e ud 
and b t h = —0.049 ± 0.003 and a t h = —0.24 ± 0.01 for e + e~ss. We see an agreement at 
the level of one standard deviation for the e + e~ case and two standard deviations for the 
case of one electron. The latter discrepancy can be a signal that our naive expectation 
for the slope is good only to a few per cent of its value. The value of the free coefficient is 
in principle sensitive to the overall normalization factor (3/4). It is however numerically 
dominated by a much larger log 10 (7 4 /7 s ) term; the actual factor of 3/4 contributes below 
10%, i.e. at the level of accuracy of the whole comparison, and is, therefore, inconclusive. 

4.2 Predictions for the 1W process 

In this subsection we present a few results of the influence of the ECS effect in KoralW on 
the liy-type process e + e~ — * e~u e ud. We use the following cut-offs (similar to Ref. ||): 

• electron angular acceptance: cos6* e > 0.997, 

• quark-antiquark invariant mass: M q q > 45 GeV, 

and the following setup: scheme; fixed W and Z widths; normal (non-screened) 
Coulomb correction; naive QCD correction; extrapolation procedure that fixes the small- 
est transfer (KeyISR=3); W branching ratios with mixing and naive QCD correction cal- 

8 The YFS exponentiation alone leads only to the exp(7/4) factor. The additional exp(7/2) comes 
from leading- logarithmic considerations. In the YFS scheme this factor is supplied order by order by 
perturbative non- infrared corrections. 



11 



culated in IBA from the CKM matrix pf; G M = 1.16639 x 1(T 5 , M z = 91.1882 GeV, 
M w = 80.419 GeV, T z = 2.4952 GeV, a s = 0.1185. 



E [GeV] 


a c+R [fb] 


1 _ oR/oC+R [%] 


1 - a c /a c+R [%] 


1 - a/a c+R [%] 


190 


87.11 ±0.22 


5.65 ±0.15 


-5.16 ±0.02 


0.81 ±0.14 


200 


103.60 ±0.26 


5.49 ±0.15 


-5.13 ±0.02 


0.67 ±0.15 


500 


807.56 ± 2.74 


4.92 ±0.21 


-4.68 ± 0.02 


0.51 ±0.21 



Table 2: The total cross-section with ECS effect (C) and running QED coupling (R). The cut-offs 
and parameter setup are as described in the text. The relative size for the ECS effect and the running 
QED coupling effect are also shown. 



In Table 2 we present the values of cross sections with the ECS effect (marked with 
C) and running of the QED coupling (marked with R) for a few energies in the LEP2 
and future LC energy range. In the respective columns we show separately the changes 
due to the ECS effect (third col.), the running of the QED coupling (fourth col.) and 
both corrections together (fifth col.). All changes are given in per cent with respect to the 
corrected cross section. One can see that the reduction of ISR due to ECS suppression 
increases the cross section for all energies, consistently with the overall effect of ISR. The 
running of the QED coupling decreases the cross section by roughly the same amount 
regardless of the energy. The relative shifts for LEP2 energies presented in Table 2 are in 
a qualitative agreement with the results of Ref . . 



4.3 Precision of the 1W cross sections 

Finally, we have to address the question of the physical precision of the above 1W cross 
sections. There are two components of the error - approximations in the ECS treatment 
and approximations related to the running QED coupling and higher order EW effects. 

• The approximations in the treatment of ECS are of the genuine non-leading type, 
i.e. a/ir with possible enhancement factors (e.g. 7r 2 ). Therefore we estimate this 
error to be below 2%. 



Following Ref. ||18|| , we estimate the approximation due to the naive treatment of 
running the QED coupling to be at most 2%, dominated by the leptonic channel 
contribution. In some specific cases, such as the semileptonic final states, this error 
can be lowered to 1% (see |L8| for details). 



The other missing EW effects we estimate at 1%, again following Ref. [18 



To summarize, in the new version 1.53 of KoralW, upon adding the above contributions 
in quadrature, we obtain an overall precison tag of the single-IT^ cross sections of 3%. 
This number can be reduced for some specific final state configurations, because of a 
smaller error contribution from the running QED coupling. The precision of the ECS 
implementation alone is of the order of 2%. Such overall theoretical precision of 3% for 
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the single-TV cross sections lies well below the expected final LEP2 experimental accuracy 
of~7%Q. 



5 Summary 

We have shown that it is possible to improve the standard ISR calculation in such a way 
that it takes into account electric charge screening and LL scale transformation for singly- 
and doubly-peripheral configuration in the production of the 4/ final states, with electron 
and/or positron in the beam pipe. This is done using a relatively simple and well-behaving 
MC weight. The method does not necessitate an explicit inclusion of the FSR. Although 
the method is not the exact implementation of the YFS exclusive exponentiation for 
the ISR+FSR, it is nevertheless closely modelled upon it. The QED running coupling 
constant for the hard process is corrected at the same time to the correct t-channel scale. 
We also present a number of numerical consistency checks and a sample prediction for the 
1W total cross section. The MC program KoralW (in its new version 1.53) in which this 
method is implemented is readily available for any interested user. It will be especially 
useful for the final LEP2 data analysis. 
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A Virtual+soft normalization weight 

In this appendix we discuss the construction of the ad hoc normalization weight of Eq. 
(|ToD and show that it ensures independence of the cross section from the dummy infrared 
cut-off e. 

Before going into details, let us stress again that we introduced the ECS ansatz by 
hand into the pure ISR-type MC algorithm in order to cure some low-angle-emission 
pathologies. Therefore it is not the purpose of this appendix to rigorously derive the 
ECS ansatz from "basic principles". Such a rigorous introduction of the ECS requires 
an entire reformulation of the MC algorithm, from ISR-type to ISR®FSR-type, based 
on the complete YFS theory and has already been done in other MC programs, such as 
BHWIDE or KKMC [0,[ll[. 

We will use the Mellin transform representation of the MC master formula as given in 



Appendix A of Ref. [|T7]]. In order to establish the notation, we will, following Ref. |p"Tf , 
briefly recall the relation of this representation to the standard Monte Carlo form of 
refs. [|T2"| , rj, 7J. Next we will show how the real emission weight is added into the master 
formula and what the matching virtual+soft compensation weight must look like to cancel 
the fictitious e dependence. 
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We start from the Mellin transform representation of the master formula as given in 
Eq. (Al) of Ref. adapted to the four-body final state (for details we refer to Ref. |T7|| ) 



a 



d 4 x 



nd 3 qi 
i=i Hi 



i x 



fpi+P2-^?i) + D 

V i=1 / 



exp (2a5 + 2aB) 



d 3 k 



ka( 3 )fJR- TZ 



1 

+ IT, 



k° 



d 3 ki d 3 k 2 d 3 k 3 



d 3 ki d 3 k 2 



-ixki—ixk 2 a{3) ( TZ TZ 



2! J kl k% 



-ixki—ixk2—ixk3r>(3)(TZ TZ 



o| / uO uO UO 
o. j \ 2 3 



(14) 



where 



d 3 k ~ 

2aB = j —S l2 {k)9{^s - k°) 



(15) 



Now the dummy IR cut-off e is introduced. Keeping in mind that the d A x integral provides 
even stronger cut-off than 9{\fs — k°) and that in the small-e limit e~ txk — > 1, one can 
rearrange the exponents 



D + 2aB = D' + B(e 

d 3 k ~ 



Be) 



k0 S 12 (k)9(e^--k ), 



(16) 



D' = j ^S 12 (k)e- ixk 9(V-s-k°)e(k -e^-). 



After expanding the D' integral and performing d 4 x integration one obtains the master 
formula in the familiar Monte Carlo form Jjj , T3j,[7|: 



a 



tUti'f (n^(*.K«-^))^(^+i*-E*-i:*.) 

n=0 ' J i=l Hl \i=l 1 J V i=l i=l 7 

exp (2aB + j ^S 12 (k)9^ - A; )) \$\pf, qf) + £ 



i>j 



Si 2 {ki)Si 2 (kj) £-^ >l S 12 (ki)S 12 {kj)S 12 (ki) 



(17) 
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Reversing the above procedure we see that the ad hoc introduction of the real emission 
weight Weqs of Eq. (||) into Eq. ( |17D amounts to the modification of the D' function: 

D'^D' R = J ^S 12 (k)w R (k)e- ixk e(yG - ky(k - e^). (18) 

Our goal is now to find the corresponding virtual+soft normalization weight that would 
exactly compensate for the change of normalization caused in the master formula by the 
replacement D' — > D' R . It is evident that such a weight should have the form 



= expp' - D' R ) = exp (J ^S 12 (k)(l - w R {k))e^ k e(k, - e^- 



(19) 



Adding and subtracting the function AU(e) of Eq. (|10D we obtain (again in the small e 
limit) 



W^ = exp(AD(e)+AU(e)) 
AD(e) = D — D R , 
d 3 k ~ 



Dr — J ^S 12 {k)w R {k)(e-^ k -e{^-s-k°) 



(20) 



Let us now summarize the situation. We have constructed the exact compensating weight 
^ exact ■ ft "would exactly restore the normalization of the master formula changed by 
the real emission weight W^ l s . At the same time, however, this weight would ruin 
the entire MC algorithm. Because of its x-dependence through the AD(e) function, the 
series of integrals over d 4 x leading to four-momenta-conserving delta functions could not 
be performed now! Therefore we have to modify the compensating weight. We have 
to do it in such a way that: (1) the master formula remains independent of dummy e 
parameter and (2) the normalization of the master formula remains as close as possible 
to the original normalization of Eq. (14). To fulfil these requirements we observe that 



the whole dependence on e is contained in the AU(e) integral. We will not modify it 
and the condition (1) will be fulfilled. The trouble-making AD(e) function is up to the 
small term (1 — w R (k)) identical to the original D function. As is known fT5| , the role 
of the D function is to compensate for the four-momentum non-conservation in the case 
of emission of multiple high energy photons. Numerically it is a very small correction, 
of the 0(a 2 ) LL type. Therefore within our accuracy we can drop this whole term and 
the condition (2) will be fulfilled within an 0(^ 2 ) accuracy. Upon including the LLST 
correction factor exp ((3/4) (j t — 7 s )j, as described in Section |2.2| , we arrive at the 
ofEq. (0). 

B Analytical integration 

We show here how to calculate analytically the AU(e) function of Eq. (|I0"D. In the polar 
variables the dk° integral decouples, as S a b(k) is proportional to l/(fc ) 2 and the ratio of 



15 



S functions is independent of the scale of k°: 

—S ab (k)(l-w R (k))=- J — J <m k F(a, b, c, d; n») = { lR - 7) log - 



ex/3/2 eV3/2 



Ir-1 = J dtt k F(a, b, c, d; tt k 
F(a, b,c,d; n k ) = (k°) 2 S ab (k) 



0,2 j- /,.\Sao( k ) + S bD(k) + S aD (k) + S bC (k) 



S ab (k) + S CD (k) 

(21) 

The evaluation of (7^ — 7) simplifies if we note that this integral is Lorentz-invariant. 
It is easiest to show for the initial form of A£/(e) given in Eq. flTUp. The only appar- 
ently Lorentz-variant part of Eq. fllPf) are the integration limits of k°. Upon a Lorentz 
transformation with an arbitrary parameter f3, the ko transforms as 

, k -(3k 1 - cos Z(fffc) 
h^k = = k A(tt k ); A(ll k ) = (22) 

1-/32 ^_^2 



and the integral AU(e) becomes 

AU(e)= J —F(a',b' } c' } d';n k ,). (23) 

A~ 1 e v ^/2 

By the rescaling transformation k' — > Ak', due to the identity Q k > — &>Ak' and the scale 
invariance of the function F(a, b, c, d; fl^), the Lorentz invariance of A/7(e) becomes trans- 
parent. 

After substituting the definition of S from Eq. (^) into Eq. ([H]), we get 

dQ k F(Q k ) = h + I 2 (24) 

Ji = / dn k (k^s ab (k) Sac{k) + SaD{k) (25) 

a 2(a6) (ac)(n k d) - (ad){n k c) 

dcosUdcp—r- -- — — : — - — — (26) 

4^ {n k a) (ab)(n k c)(n k d) + (cd){n k a)(n k b) 

dn k (krS ab (k) f^ + ^\ = h(a ~b,c~ d), (27) 
Sab{k) + S C D{k) 

where n k = k^/k°. The integral contains scalar products of k with four different four- 
momenta a, . . . ,d. With the parametrization of (n k u) = u° — u z cos 6 — u x sin # cos <p — 
u y sin 9 sin the c?0 integral in 1\ is in principle solvable. However, in the LAB frame, 
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i.e. in the CMS frame of (a, b), the (n k c) (n k d) product in the denominator would lead to a 
fourth-order polynomial and its zeros would have to be found numerically. This problem 
can be avoided by changing the Lorentz frame to the CMS(a, c), i.e. the Breit frame of 
(a, c)! The I 2 integral would similarly be evaluated in CMS(&, d) and then, due to the 
Lorentz invariance of the integrals Jj we can simply add ii(ai, . . . , d\) + I 2 (a 2 , ■ ■ ■ , d%). 
One must only remember that in both integrals the a, . . . ,d vectors became transformed 
to different values (different frames), which we indicated above as aj, a 2 , etc. 
In the CMS(a, c) we have {n k a) = a? — acosO and (n k c) = c° + c cos 8 , so that 

7i = / rfcosfl - / \ / d0 — , ns 28 

a (afo) I" L° — L z cos 9 — L x sin 6* cos — L y sin 6 sin 
COS 2?r 2 (n fc a) 7 M° — M 2 cos - M x sin cos - Mj, sin sin ^ ' 
U L = (ac)d^ - (ad)cP (30) 
M"(0) = (ab)(n k c)d^ + {cd)(n k a)lf : (31) 



The integral over d<\> can now be solved with the help of a textbook formula (cf. e.g. |EU 
2.558) 

lo + li cos0 + h sin0 
m + nil cos + m 2 sin 

_ kmi + / 2 m 2 2tt / hm, + Z 2 m 2 \ , , , 

~~ ^ ^27^2 + / 2 2 =f l + ™ 2 -i-rw 2 m o J ) m Q >m 1 +m 2 [62) 

= 27r/o/m ; mi = m 2 = 0, (33) 
to give 

T r Q a (ab) ( L X M X + LyM y , ^ - L* cos0 - (M - M z cos, 

h = GSCOS0 — — — * + 



n(n k a) ^ MJ + M 2 ^(Mq - M, cos 0) 2 - (M 2 + M 2 ) sin 2 

(34) 

There is one technical point related to Eq. ([34]), concerning the positiveness of the quantity 
under the square root in the last term. It is guaranteed, however, by the positiveness of 
the denominator M(0) ■ n k , which is a sum of dot products of time- and light-like four- 
vectors: at its minimum (with respect to 0) the M{6) ■ n k is equal to the quantity in 
question. The remaining d cos integral we perform numerically. 
The integral I 2 is evaluated analogously to Jj. 
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Figure 1: Photon angular distribution in "Mhamha scattering" [i + — ► [i~ at y/s = 5 GeV 
and muon scattering angle of 20°. The difference between left- and right-hand side plots illustrates 
the electric charge screening effect. 



20 




21 




Figure 3: Photon radiation in backward scattering, 9 = 135°. 
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Figure 4: Photon angular distribution in "Mhamha scattering" — > ^~/x + 7 at y^i = 5 GeV 
and muon scattering angle of 20° for the case of ISR only. The difference between left- and right-hand 
side plots shows the effect of the ECS correction weight. 
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-6 -4 -2 2 4 6 

V = -log lo tan(0 7 /2) 



Figure 5: log 10 of da/dlog 10 tan(0 7 /2) with (red dots) and without (blue open squares) the ECS 
correction, arbitrary units. In boxes the values of fits are shown. 
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Figure 6: log 10 of dcr/d log lo tan(0 7 /2) with the ECS correction for ev e ud (red dots) and e + e ss 
(blue open squares), arbitrary units. 
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-6 -4 -2 2 4 6 

V = -log lo tan(0 7 /2) 



Figure 7: log 10 of da/dlog 10 tan(# 7 /2) with the ECS correction for ev e ud for two values of cuts 
for electron opening angle, arbitrary units. 



26 



0.1 



0.05 



i da c /dy 
iO &10 da/dy Q t_ 



-0.05 
-0.1 

-0.15 
-0.2 



-0.25 



fit: a + b y 

a = -0.2428 + 0.0099 
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a = -0.0931+ 0.0046 

b = -0.0208+ 0.0015 
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Figure 8: log 10 of the ratio of (1 / v) (da /d log 10 v), v = 1 — s'/s, with and without the ECS 
correction for ev e ud (red dots) and e + e~ss (blue open squares) final states. 
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